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ABSTRACT 



f 



Aims. We investigate the effects of acceleration during non-linear electron-beam relaxation in magnetized plasma in the 
case of electron transport in solar flares. 

Methods. The evolution of electron distribution functions is computed using a three-dimensional particle-in-cell elec- 
tromagnetic code. Analytical estimations under simplified assumptions are made to provide comparisons. 
Results. We show that, during the non-linear evolution of the beam-plasma system, the accelerated electron population 
appears. We found that, although the electron beam loses its energy efficiently to the thermal plasma, a noticeable part 
of the electron population is accelerated. For model cases with initially monoenergetic beams in uniform plasma, we 
found that the amount of energy in the accelerated electrons above the injected beam-electron energy varies depending 
the plasma conditions and could be around 10-30% of the initial beam energy. 

Conclusions. This type of acceleration could be important for the interpretation of non-thermal electron populations 
in solar flares. Its neglect could lead to the over-estimation of accelerated electron numbers. The results emphasize 
that collective plasma effects should not be treated simply as an additional energy-loss mechanism, when hard X-ray 
emission in solar flares is interpreted, notably in the case of RHESSI data. 
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1. Introduction 

Solar flare X-ray observations provide often unique insights 
into the processes of electron acceleration and transport. 
Recent observat ions of solar flares, notably with RHESSI 
(|Lin et al.l2002|) have emphasized the high efficiency of elec- 
tron acceleration in solar flares [for a re cent review of elec- 
tron properties, see iKontar et al.l (|2011l ) and for the corre- 
sponding impli catio ns for particle transp ort in solar flares, 
iKarlickvl (l!997h and lHolman et al.1 (l2011ri ]. 

The high efficiency of electron acceleration results in 
high electron fluxes or concentrations of deka-keV elec- 
trons in solar flares, and the subsequent importance of 
collective effects to particle transport in the solar atmo- 
sphere. The presence of large number of energetic elec- 
trons in coronal loops could trigger a number of instabil- 
ities and generate plasma waves, which in turn affect the 
transport of energetic particles from the acceleration re- 
gion down to the chromosphere. It has been shown that 
accounting for these collective effects could affect the in- 
terpretation of hard X-ray spectra and lead to additional 
observational consequences that are essential to the study 
of solar flares. For example, the inclusion of Langmuir 
wave generation in the treatment of spatially localized 
electron beams (Hannah ct al. 2009) prevents the forma- 
tion of a pronounced low-ener gy cut-off th a t appears in 
purel y collisional models ( e.g. Brown 19711 : iBrown et al.l 
l2002tl . Weibel instability (We ibell Il959h can quickly in- 



crease the velocities in the direction perpendicular to the 
beam propagation and affect the observed X-ray aniso tropy 
(jKarlickv fc Bartall2009HKarlickv fc Kasparovall20090 . The 
presence of Langmuir waves in a flaring loop could also re- 
sult in plasma emissio n (e.g.lVlahos fc Papadopoulosl ll979l: 
lEmslie fe Smith|[l98 4; Ham ilton fc Petrosianll987[ ). provid- 
ing additional constraints on non-thermal electron popula- 
tions in solar flares. 

In laboratory plasma experiments, collision-less effects 
involving various instabilities have been demonstrated to 
play a key role in electron transport. In experiments 
to study the beam-plasma interaction, the appearance 
of electrons with energies exceeding that o f the injected 
beam energy was not e d in early studies (e.g. iBerezin et al.l 
11964 iFainberd [1961 iKovalenkol fl983l ). These above-the- 
injected energy electrons are normally conn ected to either 
the presence of plasma inhomogeneities (IVedenov et al. 
1967t iRvutovl 119691: iNishikawa fc Rvutovl 119761: lEscande 
1979T ) or the nonlinear effects of wave-particle inter- 



actions (|Breizman fc Rvutovl 119751 : iTimofeev fc Terekhovl 
I2010T ) . Furthermore, the influence of plasma inhomo- 
geneities or wave-wave interactions on beam-plasma in- 
stability in the solar context has been studied exten- 
sively in conne ction to the theory of solar radio- type 



III bursts (e.g. Muschiett et al. 19851: Tvshetskiv et al 



20071 : IKontar fc Reidll2009l: lTsiklauril l2010l: iReid fc Kontarl 
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Table 1. Model parameters and the fraction (FR) of the 
beam energy in the electrons with energies greater than the 
initial-beam electron energy. 



Model 


mi/m c 


n b /n c 


Vb/C 




FR (%) 


A 


16 


1/8 


0.666 


0.0 


10 


B 


16 


1/8 


0.666 


0.1 


10 


C 


16 


1/8 


0.666 


0.5 


22 


D 


16 


1/8 


0.666 


0.7 


28 


E 


16 


1/8 


0.666 


1.0 


29 


F 


16 


1/8 


0.666 


1.3 


27 


G 


16 


1/40 


0.666 


0.0 


12 


H 


16 


1/40 


0.666 


1.0 


28 



sands of seconds. Therefore, modeling efforts in solar flare 
physics have focused either on relatively long observational 
timescales by either ignoring short-timescales or addressing 
the micro- physics on s cales smaller than the observational 



e.g. lBredl2009t iBingham et alj[200Tl: lEliasson et all 



scale s 

20061: iHewett et all 119911 _ . 

Rose et al.ll2002t ISakai et al.ll2006l ISaito k Sakaill2004 
Kontar et al.l (|2012l ) showed that the £;-spectrum evo- 



Lee et all 120081 TlVlessmerl 12002 



lution of beam-generated Langmuir waves in collisional 
plasma could result in substantial energy gain by high en- 
ergy electrons or the effective acceleration of electrons of the 
same beam above 20 keV for solar flare conditions. Thus, 
if collisional relaxation is assumed, the number of ener- 
getic electrons inferred from X-ray spectra could be overes- 
timated and may lead to an apparently large number of ac- 
celera ted electrons. However, the treatment of lKontar et all 
(I20l1 is a one-dimensional analysis based on weak turbu- 
lence theory. They showed that a positive density gradient, 
externally excited density fluctuations, and the three-wave 
interaction involving ion-sound mode could lead to an ef- 
fective acceleration of beam electrons. However, the role 
of the guiding magnetic field as well as three-dimensional 
(3D) aspects of beam-plasma interaction have not yet been 
addressed. 

This paper investigates the non-linear evolution of non- 
thermal electron beams in a plasm a. Using t he pa rticle- 
in-cell (PIC) model develop ed in iKarlickvi (|2009h and 
iKarlickv'fe Kasp aroval (|2009l ). we investigate the effects of 
the acceleration of electrons in the beam during a non-linear 
stage of beam-plasma instability for various values of the 
guiding magnetic field. We show that about ~ 10 — 30% of 
electrons are accelerated to the energies greater than the 
energies at which they were injected. The presence of a 
guiding field increases the number of electrons accelerated 
during a beam-plasma interaction. The results show the 
appearance of accelerated electrons and highlight that col- 
lective effects can lead to not only additional energy losses 
for lower energy electrons but the acceleration of some elec- 
trons above the initial energy of electrons. As the hard X- 
ray spectra in solar flares is normally due to deka-keV elec- 
trons, this acceleration effect should be taken into account 
when the number of accelerated electrons is estimated. In 
this paper, we compared these results with analytical esti- 
mates. 



2. Simulation model 

It is commonly accepted that electrons in solar flares are ac- 
celerated at low corona heights by primary energy-release 
processes. They propagate along the magnetic field lines 



as electron beams downwards to loop footpoints, where 
they generate hard X-ray emission. The rapid variations (of 
timescales ~ 45 msec) in this hard X-ray emission observed 
in some events indicate that the acceleration as well as the 
electron beam flux can may themselves rapidly varying (e.g. 
iKiplinger et~aTlll983l:lAschwanden et alj fl995). This means 
that the effect of fast propagation needs to be considered in 
the evolution of these electron beams. Assuming that the 
electrons have a power-law distribution in the acceleration 
region. Owing to the propagation, fast electrons could then 
overtake slower ones at some distance from the acceleration 
site (e.g. lHannah et aL1l2009h . forming an unstable distri- 
bution that can be approximated by the mono-energetic 
beam. 

For our study, we used a 3D (three spatial and three 
velocity compo nents) relativistic electromagnetic PIC code 
(|Karlickvll2009D . The system sizes are L x — 45A, L y = 45A, 
and L z — 600A, where A is the grid size. 

We initiated a spatially homogeneous electron-proton 
plasma with the proton-electron mass ratio m p /m e =16 
(models A-H in Table 1). This ratio was chosen to 
shorten the computational times and keep the proton skin- 
depth shorter than the dimensions of the numerical box. 
Nevertheless, the ratio was still sufficient to clearly sepa- 
rate the dynamics of electrons and protons. The electron 
thermal velocity was vt c — 0.06 c, where c is the speed of 
light. In all models, 160 electrons and 160 protons per cube 
grid were used. The electron plasma frequency was w pc = 
0.05 (At) -1 (At = 1 is the time step) and the electron 
Debye length was Ad = VTc/^Jpc = 0.6 A. The electron and 
proton skin-depths were A cc = c/io pc — 10 A and A cp — 
c/ojpp = 40 A (where w pp is the proton plasma frequency), 
respectively. 

We then included a mono-energetic beam that was ho- 
mogeneous throughout the numerical box. To keep the to- 
tal current zero in these models in their initial states, we 
introduced an initial return current by shifting the back- 
ground plasma electrons in the velocity space according to 
the relation = —v^n^/rie, where v\, is the velocity of the 
electron beam, and rib and n c are the beam and background 
plasma densities; for an exam ple of this type of initializa- 
tion, see lNiemiec et all (|2008f ). In principle, it is possible to 
start from a zero initial return current. However, owing to 
inductive effects included in the used 3-D electromagnetic 
code, the current starts to oscillate at the electron plasma 
frequency with the amplitude of the stabilizing return cur- 
rent. Therefore in these simulations, we recom mend start- 
ing f rom an already stabilized return current (Nishikawa 
2008). Although these current oscillations can influence the 
studied acceleration, for simplification, we did not take this 
effect into account. We similarly neglected any effects of the 
return-current formation (electrostatic ones) on the front of 
the propagating beams. 

The beam velocity was chosen to be v-^/c = 0.666 in the 
z direction. The ratio of the beam to the plasma densities 
was chosen to show the details of the acceleration process 
rih/n e = 1/8, which is a rather high value (however still 
realistic). The lower values of nb/n c = 1/40 , see Table 1, 
for comparison, were used in models G and H. 

To investigate the influence of the magnetic field in the 
models, we consider several values of the background mag- 
netic field, so that the ratio of the electron-cyclotron to 
electron-plasma frequencies (w cc /cj po ) is 0.0, 0.1, 0.5, 0.7, 
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Fig. 1. The electron velocity distributions for model E at four different times: at the initial state (a), at u) pc t = 40 (b), 
at LO pe t = 140 (c), and oj pc t — 200 (d). Crosses correspond to f(v z ), dotted and dashed lines display f(v x ) and f(v y ), 
respectively. Note that f{v x ) and f(v y ) overlap. The vertical line in the part a) at v/c = 0.666 denotes the monoenergetic 
electron beam. 



1.0, and 1.3 (see Table 1). In all models, the periodic bound- 
ary conditions were used. 

3. Results of numerical simulations 

Using the above described model, we performed eight runs 
(A-H) using the parameters given in Table 1. As an example 
(see Fig. [T]) , a time evolution of the electron velocity distri- 
bution for model E is shown at four different times: the ini- 
tial state (a), uj pc t = 40 (b), uj pc t — 140 (c), and uj pc t — 200 
(d). Crosses correspond to f(v z ), and dotted and dashed 
lines display f(v x ) and /(%), respectively. The vertical 
line in the initial state at v/c = 0.666 denotes the mono- 
energetic electron beam. As can be seen here, a plateau 
forms in the velocity space between the velocity of the ini- 
tial beam and that of the background plasma electrons. 
It is produced by the so-called quasi-linear relaxation, in 
which the electron beam generates the Langmuir waves (see 
Fig. [5]). The Langmuir waves were initially generated at 



the fc-wave vector corresponding to the resonance condition 
for the beam-plasma instability, i.e. at lj = k z Vb (see the 
peak in the Langmuir wave energy-spectrum in Fig. [2^), 
where uj is the frequency of t he Langmuir waves and 
is the electron beam velocity (jMikhailovskiil [1975F) . These 
Langmuir waves then decay and merge mainly in accor- 
dance with the th r ee wa ve interactions (for details, see e.g. 
iBarta fc Karlicky] (|2000h ') and evolve in the k- vector space 
(see the enhanced Langmuir wave spectrum in Fig.[2j3,c,d). 
The Langmuir waves, on the other hand, scatter the beam 
electrons and heat the background plasma. During this pro- 
cess, the plateau in the electron distribution function f(v z ) 
is formed; (see detailed discussion of the quasi-linear relax- 
ation in the following section). However, during this pro- 
cess some beam electrons, owing to an interaction with 
the Langmuir waves, obtain an energy that is higher than 
the electron energy in the initial beam. This is shown in 
Fig. O where a time evolution of the electron energy dis- 
tributions for model E at four different times is presented: 
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Fig. 2. The Langmuir wave energy for model E in the fc-space at four different times: uj pc t = 40 (a), uj pc t = 60 (b), ui pe t 
= f40 (c), and uj pc t = 200 (d) (solid lines). For comparison in each panel, the initial Langmuir wave energy is added 
(dashed line). 



at (a) uj pc t = 40 (b) uj pc t = 60 , (c) uj pc t = 140, and (d) 
uj pc t = 200. For comparison, we show in each panel in this 
figure the initial electron-plasma distribution together with 
the initial mono-energetic beam (dashed lines). Here, it can 
clearly be seen that some electrons have the energies higher 
than those of the initial beam. 

Although in the initial state, we started from a uniform 
background plasma, very soon (owing to the dense beam) 
strong density fluctuations appeared (see Fig. [5]). During 
the evolution, their characteristic lengths became longer 
(compare Figs. [5^, andlSJa). In the early stages of evolution, 
the associated electric field densities are not correlated with 
the density fluctuations. On the other hand, in the later 
stages of evolution, e.g., at uj pc t = 190 it can be seen that 
the electric fields (Langmuir waves) start to be trapped in 
density depressions (see e.g. the electric field densities at z, 
equal to 300 A and 430 A respectively). 

To investigate the energetics of the accelerated elec- 
trons, we computed their energy (expressed in fraction of 



the initial beam energy) above some selected energy lev- 
els. The time evolution of these fractions for the energies 
higher than E/mc 2 — 0.25 and E/mc 2 — 0.3, for model E 
are shown in Fig. [5J The initial energy of the beam elec- 
tron was E/mc 2 = 0.22178. In this figure, for comparison 
the time evolution of the maximum energy of the electron 
(E/mc 2 ) is plotted. The fractions firstly increase in time, 
owing to the some time taken for acceleration of these elec- 
trons to the selected energy levels. The fractions then slowly 
decrease to some saturation level, in agreement with the 
saturation of the quasi-linear relaxation process. The frac- 
tions generally decrease with the increase in the energy in- 
terval between the beam energy and selected energy level. 
We also computed the fractions of these accelerated elec- 
trons for the energy levels just above the beam energy. For 
all computational models (A-H) at the time uj pc t = 200, 
these fractions are summarized as percentages in the last 
column of Table 1 . These fractions increase as the magnetic 
field increase, e.g. from FR = 10 % for ui cc /uj pc — (model 
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Fig. 3. The electron energy distributions for model E at four different times: uj pc t — 40 (a), w pc t = 60 (b), oj pc t = 140 
(c), and w pe i = 200 (d) (solid lines). For comparison in each panel, the initial electron-plasma distribution together with 
the initial monoenergetic beam are added (dashed lines). 



A) to FR = 29 % for w cc /u pc = 1 (model E). For even 
greater values of uj cc /lo pc , the fractions start to decrease 
(see FR = 27 % for ui cc /u pc = 1.3 in model F). Similar re- 
sults were found for models G and H, where we considered 
a beam density lower than in models A-F (see Table 1). 
As sh own a nd analyzed in the papers of | Karlickv fe Bartal 
(|2009f) and iKarlickv fc Kasparova! (|2009t ) for the models 
with weak magnetic fields (models A-B), the Weibel in- 
stability is significant and strongly influences the result- 
ing electron distribution function. On the other hand, the 
strong magnetic field (e.g. model E) reduces the role of the 
Weibel instability (see also the discussion and Fig. [5]). 



4. Quasi-linear relaxation - analytical estimates 

While the analytical treatment of a 3D non-linear beam- 
plasma system is impossible, the efficiency of accelera- 
tion can be estimated in the quasi-linear limit. In this 
limit, the evolution of energeti c electrons is described by 
two coupled k i netic equations ( Vedenov fc Velikhovi ll963: 
iVedenov et all [19671 lTsvtovichlIT995r For a weakly mag- 
netized plasma, when the field is strong enough to provide 
one-dimensional electron dynamics, the electron beam evo- 
lution can be described using standard quasi-linear theory 



dt 



v dv 



As an illustration, we present in Fig. [4] the electron en- 
ergy distributions at uj pc t = 200 for all runs, i.e. for uj cc /uj pc dWk du> pe {x) dWu 



0.0, 0.1, 0.5, 0.7, 1.0, and 1.3, respectively (models A-F, 
Tabic 1). The number of electrons accelerated above the ini- 
tial electron-beam energy (expressed by the vertical dashed 
line) corresponds to the fractions FR in Table 1. 



dt 



dx dk 



TUJp e ^ e Of 

n e k 2 dv 



(1) 



(2) 



where f(v,t) is the electron distribution function, Wk is 
the spectral energy density of the Langmuir waves, u> pe 
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Fig. 4. The electron energy distributions (solid lines) at uj pc t = 200 as a function of the magnetic field in models A-F 
with uj cc /uj pc = 0.0, 0.1, 0.5, 0.7, 1.0, and 1.3, respectively. For comparison in each panel we plot the initial electron 
plasma distribution together with the initial monoenergetic beam (dashed lines). 



1 .0 

0.5 
0.0 



[ a) 









100 200 300 400 500 

z (A) 



BOO 



2.0 p 
1 .5 
1 .0 -' 




100 200 300 400 500 

z (A) 




Fig. 5. Normalized plasma density (solid line) and the elec- 
tric field density (dashed line) in the z-direction along the 
line with x = y = 22A for model E at two times: Lu pc t — 
60 (a) and u pc t — 190 (b). The electric field density is ex- 
pressed in arbitrary units, but the electric field density at 
(b) is multiplied by a factor of 50. 



Fig. 6. Fractions of the beam energy in electrons with en- 
ergies greater than E/mc 2 = 0.25 (dashed line) and greater 
than E/mc 2 = 0.3 (dotted line), for model E. The initial 
energy of the beam electron is E/mc 2 — 0.22. For compar- 
ison, we plot the time evolution of the maximum energy of 
the electron (solid line). 



is an electron plasma density, and e and m are the elec- 
tron charge and mass. The equations describe the resonance 
(uj pe = kv) interaction of the electron beam with the sur- 
rounding plasma via the generation of plasma waves and in- 



clude plasma inhom ogeneity effects (see, e.g . IVedenov et al.l 
119671: iRvutovlllflfia iKontar fc Pecselill2002li . 
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Fig. 7. The X-ray spectrum for model E at u) pc t = 200 
(solid line) and the initial state (dashed line). 



When the plasma is uniform 



dm pc (x) 
dx 



0, the station- 



ary solution for the initially unstable beam distribution 
/(v, t = 0) = go(v), v < Vh of the coupled quasi - linear 
equations is well-known (|Vedenov & Velikhovl [l963h . The 
electron distribution function has the form of a plateau 



f{v,t^r oo) = — , 

Vb 



v < Vb, 



(3) 
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and the spectral energy density of plasma waves Wk can be 
found from Equations (|1I2[) using Equation @ 

W<yk = ~V' 1 ~* ^ = Iq {y b ~ dV ' V < Vbl ^^ig- 8 - Comparison of the electron distribution functions in 

model A (without the magnetic field) and model E {uj ce /ui pe 
= 1.0). Crosses correspond to f{v z ) and dashed lines display 
overlapping f(v x ) and f(v y ). 



thus the spectral energy density becomes 



v 



mn b 4 
oo) = v , v < v b 

UJ pe Vb 



for go(v) = rib5(v — Vb) and takes the form 



v 



mni, q 
v 

UJpeVb 



V 

1 

Vb 



V < Vb 



(5) 



(G) 



for go(v) — 2nbvjv\, v < Vb, (see lKontanl2001l for details). 
The total energy density of plasma waves from Equation 
© is U L = 2/3 x mn b vl/2 and U L = 1/3 x mn b vl/2 for 
the second case given by Equation ©. 

However, when the plasma has a positive density gradi- 
> 0, which corresponds to electrons propagating 



ent 



dx 



into the region of higher density (as in the standard flare 
scenario), the Langmuir waves slowly evolve towards larger- 
phase velocities or smaller wavenumbers k. As Langmuir 
wave packets propagate within the plasma, t he total energy 
of the wave packet u>(k, x) must be constant (|Vedenov et al.l 
119671) . which requires that the wavenumber change be neg- 
ative for a positive plasma-density gradient, such that 



Ak ~ ^AJ-M. 

ox 



(7) 



Therefore, the motion of Langmuir wave over the time 
range At results in a decrease in the wavenumber from k 
to k - ^SlAt. The waves shifted to higher v (smaller 
k) can be effectively re-absorbed by the beam, which leads 
to an acceleration of the electrons and the formation of 
an extended plateau above v > Vb- Using the conservation 
of energy, and assuming that all Langmuir waves are re- 
absorbed by the beam owing to the plasma inhomogeneity, 
W(k = oj/v, t — > oo) = 0, one finds that for the energy 

/•OO 

f{v,t -> oo)v 2 dv = / f(v,t = Q)v 2 dv (8) 
Jo 

the electron distribution has electrons with velocity v > v b 

Tib 

/(u,i->oo) = , v<v max , (9) 

Umax 

where the new maximum velocity is from Equation (JSJ) 
3 



rib Jq 



go(v)v dv. 



(10) 
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Numerical solutions of equations ([T][2]) (see iKontarl I2001L 
for details) show that the value of the density gradient 
mostly affects the rate of the extended plateau formation 
and that the final state of the initially unstable distribution 
g (v) = 2n b v/v b 2 , v < v b , gives v max = y/ij2v b . In the case 
of g (v) = n b 5(v - v b ), one finds that = V3v b . 

In the case of non-linear interactions and various den- 
sity fluctuations, the energy exchange betwee n the electrons 
and pl asma waves becomes more complicated. iKontar et aLl 
(|2012t) used numerical simulations to estimate the role of 
these effects. Here, assuming that small-scale density fluc- 
tuations in the plasma generated during the beam-plasma 
instability are random with zero mean, the evolution of 
Langmuir waves in plasma can be app roximated as a sym - 
metric diffusion in /c-space (see Eq 14 in lKontar et alj20 12). 
This leads to equal numbers of plasma waves spreading 
towards smaller and larger the phase velocities. The for- 
mer waves will be Landau-absorbed by the thermal plasma, 
while the latter will be absorbed by the beam resulting in 
electron acceleration in the tail of distribution. Therefore, 
for the Langmuir wave spectrum flattened in fc-space, (see 
Figure [2]), these simplistic arguments suggest that only half 
of the Langmuir wave energy given by Equations ([SJ ^ will 
be reabsorbed back. The maximum velocity then becomes 
Vmax = V%v b and the energy of the electrons with v > v b 
becomes 

m f Vma:c n b 2 
U(v > Vb, t — >• oo) = — / v av 

^ J Vb ^max 

-0.43-^, (11) 

for go(v) = n b 5(v — v b ). We note that these estimates are 
rather close to the numbers inferred from the 3D PIC sim- 
ulations presented in Table 1. It is worth noting that the 
collisional relaxation of the electron power-l aw spectrum, 
which is initially stable was considered by IKontar et al.l 
(|2012f ). while in this paper we treat a 'classical' case of 
beam-plasma instability. 

5. Discussion and conclusions 

We have performed a number of 3D PIC simulations of 
the beam-plasma instability with monoenergetic beams and 
have shown that during relaxation a population of electrons 
with velocities exceeding those of the injected electrons ap- 
pears. The energy of these electrons is around 10 — 30% of 
the initial beam energy. 

Using PIC simulations, it is difficult to predict the 
long-term time evolution of these processes in solar flares. 
However, as shown in Fig. [7] this effect is indicated by the 
high-energy limit of the X-ray spectrum, namely that the 
accelerated high-energy electrons shift the X-ray spectrum 
to higher energies. This result together with the radio di- 
agnostics can be used for an estimation of these accelera- 
tion processes. For example, if we take the dm-spikes as a 
radio signature of t he acceleration process in solar flares 
(|Guedel et al.l Il99ll ). then the advanced theory of these 
bursts can be used to estimate the electron distribution 
function at the acceleration site. Comparing this function 
with that determined from the hard X-ray spectrum at the 
flare footpoints, the acceleration efficiency can then be es- 
timated. 



We have found that the increasing magnetic field 
strength leads to a larger fraction of accelerated electrons. 
Therefore, we decided to compare the electron distribution 
functions of the cases with and without a magnetic field 
(in models A and E, Fig. El ). As p resented and analyzed 
in iKarlickv fe; Bartal f|2009h and Karlickv & Kasparova 
(l200l 7 the main difference in both cases is caused by the 
Weibel instability. In model E, the Weibel instability is re- 
duced, while in model A (without the magnetic field) the 
Weibel instability transfers the beam energy to a heat- 
ing of mainly perpendicular components of the background 
plasma. Thus, in the case without the magnetic field, not 
only the bump-on-tail instability but also the Weibel insta- 
bility operates and less energy (than in model E) is trans- 
ferred to the Langmuir waves, leading to a weaker acceler- 
ation. 

We considered two values of the ratio of the beam to 
background plasma densities n^/n e = 1/8 and 1/40, which 
imply that the return-current electron speeds are Vd — 0.083 
c and Vd = 0.016 c (where c is the speed of light), re- 
spectively. This corresponds to two regimes of the return- 
current electron speed either greater or lower than the ther- 
mal plasma velocity, which is in our model vt c — 0.06 c. 
In both cases, the percentage of the energy in accelerated 
electrons is similar. The case where rib/n e = 1/8, i.e. the 
case with the drift speed greater than the initial thermal 
velocity, should be unstable for the Buneman instability. 
However, the time of evolution considered in the present 
study is shorter than the time required to develop such an 
instability. For much longer term evolutions of similar sys- 
tems in this regime computed with a 1-D Vlasov code, the 
formation of weak double layers was found and proposed as 
an expla nation of broken power-law X-ray spectra in solar 
flares bv lLee et aD (|2008H (see also Karlicky 2012). 

We note that the size of our numerical box 
(45 A x 45 A x 600 A) is limited by the memory and speed lim- 
its of our computer. The limited number of spatial grids lim- 
its a number of grids in the fc-vector space. This limitation 
can influence the wave- wave and wave-particle interactions 
in this acceleration process. This is especially important in 
the linear regime of these interactions, which would corre- 
spond here to the cases with low density beams. However, in 
the present study we considered very dense beams (n b /n e 
= 1/8 and 1/40), which generate strong density fluctua- 
tions ensuring that these interactions are far from linear 
regimes (We note that in the interpretation of the hard X- 
ray emission the dense electron beams are often required). 
We made additional tests for model E. We used the same 
parameters as in model E, but we varied the size of the 
numerical system. We obtained similar results. For the sys- 
tem size (20Ax20Axl500A) at the time oj pc t = 200, the 
fraction FR is 26 % and for the size (lOAx 10Ax3000A) at 
the same time the fraction FR = 25 %; compare with the 
fraction FR = 29 % for model E in Table 1. Nevertheless, 
to make these results more precise we plan to repeat these 
computations in larger numerical boxes on a more powerful 
computer. 

Comparing the numerical and analytical treatment, we 
propose that for this type of acceleration the density fluctu- 
ations and non-linear wave- wave interactions are essential. 
While in a strictly uniform plasma this acceleration is im- 
possible (see the analytical estimations), in real conditions 
with sufficiently dense electron beams some beam electrons 
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are accelerated to energies greater than the initial ones via 
non-linear wave interactions and density fluctuations. 

In solar flares, high non-thermal electron fluxes are of- 
ten required to explain the observed X-ray emission, which 
provide the suitable conditions for fast beam-plasma insta- 
bility. As we have shown, the beam-plasma instability gen- 
erates the Langmuir wave turbulence. The acceleration of 
electrons occurs owing to the /c-space evolution of Langmuir 
waves, in which not only the Langmuir waves with phase 
velocities smaller than the initial beam velocity are gener- 
ated (as in the case of beam relaxation in uniform plasma) 
but also those with higher phase velocities. The process 
of electron acceleration is fast, and occurs on a timescale 
much shorter than the electron transport time from an ac- 
celeration site to the dense chromospheric region. We there- 
fore expect that during these processes the tail of the dis- 
tribution extends towards higher energies. The number of 
high-energy electrons above some energy level increases and 
thus influences the hard X-ray spectrum. However, the to- 
tal number of beam electrons is conserved and the process 
only redistributes the beam energy. These results empha- 
size that the transport of electrons should not be treated 
using a single particle description, and that the collective 
effects produce not only energy loss, but an effective accel- 
eration of electrons. The analysis of a hard X-ray spectrum 
assuming only collisional losses could therefore infer that a 
larger number of electrons than actually is initially acceler- 
ated owing to the additional in-flight acceleration. 
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